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Abstract: Aiming at the observation of cosmic-ray chemical composition at the “knee” energy region, we have been 
developing a new type air-shower core detector (YAC, Yangbajing Air shower Core detector array) to be set up at 
Yangbajing (90.522° E, 30.102° N, 4300 m above sea level, atmospheric depth: 606 g/m 2 ) in Tibet, China. YAC 
works together with the Tibet air-shower array (Tibet-Ill) and an underground water cherenkov muon detector array 
(MD) as a hybrid experiment. Each YAC detector unit consists of lead plates of 3.5 cm thick and a scintillation 
counter which detects the burst size induced by high energy particles in the air-shower cores. The burst size can 
be measured from 1 MIP (Minimum Ionization Particle) to 10® MIPs. The first phase of this experiment, named 
“YAC-I”, consists of 16 YAC detectors each having the size 40 cm x 50 cm and distributing in a grid with an 
effective area of 10 m 2 . YAC-I is used to check hadronic interaction models. The second phase of the experiment, 
called “YAC-II”, consists of 124 YAC detectors with coverage about 500 m 2 . The inner 100 detectors of 80 cm x 
50 cm each are deployed in a 10 x 10 matrix from with a 1.9 m separation and the outer 24 detectors of 100 cm x 
50 cm each are distributed around them to reject non-core events whose shower cores are far from the YAC-II array. 
YAC-II is used to study the primary cosmic-ray composition, in particular, to obtain the energy spectra of proton, 
helium and iron nuclei between 5xl0 13 eV and 10 16 eV covering the “knee” and also being connected with direct 
observations at energies around 100 TeV. We present the design and performance of YAC-II in this paper. 
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1 Introduction 

The all-particle energy spectrum of primary cosmic 
rays is well expressed by a power law dN/dE oc E 1 over 
many orders of magnitude, with 7 changing sharply from 
2.7 to 3.1 at about 4 PeV [l.] . The break of the all-particle 
energy spectrum is called the “knee” [ 3 , S], and the corre¬ 
sponding energy range (10 15 eV - 10 16 eV) is called “knee 
region”. Although the existence of the knee is well con¬ 


firmed by many experiments, there still exist debates on 
its origin. In order to clarify the sharpness of the knee, 
precise measurements as possible of the primary spectra 
of individual components including the knee will be es¬ 
sentially important. The best ways to study the chemical 
compositions are direct measurements of primary cosmic 
rays on board balloons or satellites, but its energy range 
with sufficient statistics is limited to 10 14 eV because of 
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limited exposure time and small detective area. So the 
task of studying chemical components of knee region still 
relies on the ground-based indirect measurements. 

The study of the mass composition of primary cosmic 
rays around the knee has been made in the last decade 
by the Tibet-EC (emulsion chamber) experiment using 
high threshold (> TeV) air-shower core detector sensi¬ 
tive to the primary particle mass Q. It was aimed to 
separate air-shower events induced by primary light com¬ 
ponent such as protons and helium nuclei. The observa¬ 
tion shows low intensities of proton and helium spectra 
which amount to less than 30% of all-particle spectrum, 
suggesting that the knee is dominated by nuclei heav¬ 
ier than helium. Demerits of this experiment are: l)tlre 
statistics are limited due to the selection of high-energy 
core events; 2)the reconstruction of the primary energy 
spectrum is based on the AS simulations in which the 
hadronic interaction model is not fully established yet, 
although the model dependence of AS core events at high 
observation level is at most few tens percent as already 
reported Q. To overcome these problems, we have re¬ 
cently made the upgrade of the Tibet-EC experiment and 
started a new low threshold core detector named YAC 
(Yangbajing Air shower Core detector), which has a ca¬ 
pability of observing the core events with high statistics 
as well as testing the interaction models. 

In order to further expand the energy region of the 
Tibet-EC experiment, YAC detector has been developed 
to lower the threshold energy as possible with wider dy¬ 
namic range. The detection threshold energy of the YAC 
is able to set at several 10 times lower than EC (about 
300 GeV corresponding to the primary energy of several 
times 10 TeV) by adopting the scintillator instead of the 
X-ray film for the detection of cascade showers induced 
in the lead plate by high energy AS core particles, and 
the wide dynamic range of 1-10 6 MIPs (Minimum Ion¬ 
ization Particles) for the burst size detection is realized 
by installing 2 PMTs (high-gain and low-gain PMTs). 
Using the wave length shifting fiber to collect the scin¬ 
tillating light improves the geometrical uniformity. This 
new experimental condition improves the statistics of the 
high energy core events by a factor of 100, when com¬ 
pared with the Tibet-EC experiment . The new hybrid 
experiment aims to observe the energy spectra of proton, 
helium and iron whose energy range will overlap with di¬ 
rect observations at lower energies such as CREAM Si 
and ATIC Q, and Tibet-EC experiment at higher ener¬ 
gies. Furthermore, we may add the underground muon 
detector (MD) to this experiment to increase the mass 
separation power of primary particles. Hence, the new 
Tibet hybrid experiment consists of YAC array, Tibet-Ill 
array [l| and the underground MD array (Fig. CO- 

The design of the YAC detector is described in sec¬ 
tion 2. YAC experiment is scheduled in two steps, named 


“YAC-1” and “YAC-II”. YAC-I is a small array con¬ 
sisting of 16 prototype detectors covering area of about 

10 m 2 located near the center of the Tibet-Ill array. 
This array will operate for a few months to observe 
the AS core events of primary energies around 10 14 eV 
where the mass composition of primary cosmic rays is 
fairly known by direct observations |5|-|7|. Therefore, 
the role of YAC-I is to test the interaction models cur¬ 
rently used by Monte Carlo (MC) simulations such as 
QGSJET01 (i, SIBYLL2.1 @, DPMJET 0 and lately 
EPOS-LHC [Hi, QGSJETII-04 0. YAC-II is an ar¬ 
ray of 124 detectors covering area of ~500 m 2 to obtain 
proton, helium and iron spectra with high statistics in 
energy range between 5xl0 13 eV and 10 16 eV, which will 
be smoothly connected to those by directs observations 
in lower energy region. In this paper, the design of YAC- 

11 and its performance are described. 

2 YAC-II Experiment 

2.1 The design of YAC-II array 

The YAC array is set up at Yangbajing (90.522° E, 
30.102° N, 4300 m above sea level, atmospheric depth: 
606 g/m 2 ) in Tibet, China. This array works together 
with the Tibet-Ill and underground MD array as a hy¬ 
brid experiment, as shown in Fig. [T] 



x(m) 


Figure 1. Schematic view of (YAC-II+Tibet- 
III+MD) array. The YAC-II array contains 124 
core detectors (blue open squares) and is located 
at near the center of Tibet-Ill (~50000 m 2 ). The 
Tibet-Ill array consists of 576 fast timing (FT) 
detectors (black open squares), 28 density (D) 
detectors (red half open squares) and 185 FT/D 
detectors (red solid squares) at periphery. The 
MD array consists of 12 underground water pools 
(gray solid squares), each of which has 16 units. 5 
water pools (~4500 m 2 ) in blue wire frame have 
been built. 
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YAC-II consists of 124 core detectors covering the 
area of about 500 m 2 , as shown in Fig. [2] The inner 100 
detectors 80 cm x 50 cm each are deployed in a 10 x 10 
matrix from with a 1.9 m separation and the outer 24 
detectors 100 cm x 50 cm each are distributed around 
them to reject non-core events whose shower cores are 
far from the YAC-II array. The main purpose of YAC-II 
is to observe the energy spectra of light component as 
well as iron nuclei in the primary cosmic rays between 
5xl0 13 eV and 10 16 eV as mentioned in Section 1 


YAC-II array 



Figure 2. Left: The magnified schematic view of 
the YAC-II array (~500 m 2 ). The detectors of 
this array are classified to 2 types: the inner 100 
detectors of 0.4 m 2 each (blue solid squares) and 
the outer 24 detectors of 0.5 m 2 each (gray solid 
squares). Right: The schematic view of each de¬ 
tector comprising a lead plate of 3.5 cm thickness, 
a iron plate of 0.9 cm thickness and a plastic scin¬ 
tillator of 1 cm thickness. 



Figure 3. The view of the scintillation detector. 
Top: the top view and Bottom: the back view. 
The scintillation detector unit consists of 20 pieces 
scintillator whose size is 50 cm x 4 cm x 1 cm 
and light-isolated from each other using reflecting 
material. All the fibers are adjusted to the same 
length of 110 cm. Both of 20 fibers attached to the 
low-gain PMT (R5325) and 40 fibers attached to 
high-gain PMT (R4125) are installed in the center 
of the scintillator with equal distance. 


2.2 The design of YAC-II detector 


The design of YAC-II detector essentially followed 
that described in 13] but in this work some improve¬ 
ments were made as follows: 

(1) Lead plates: A lead layer is used to select high 
energy particles in the AS core in the energy range from 
several GeV to several 10 TeV, as shown in Fig. [2] It 
is found from an optimization calculation that the lead 
layer of thickness 3.5 cm (6.3 radiation length) meets 
the measurements in both lower and higher energies. A 
0.9 cm thick iron plate placed in the middle of scintil¬ 
lator and lead plates is used to support the weight of 
lead plates, as shown in Fig. [2] High energy electro¬ 
magnetic particles near the AS axis develop into cascade 
showers in the lead absorber and these shower particles 
enter to the scintillator counter. Here, the number of 
shower particles detected in the scintillator is defined as 
the burst size (N b ). The scintillation light produced by 
shower electrons in the scintillator below the lead layer is 
transmitted to PMTs via the wave length shifting fibers 
(WLSF BCF-92, SAINT-GOBAIN, round cross section 
with 1.5 mm diameter), as shown in Fig. [3] The deter¬ 
mination of the burst size ( N b ) is calibrated by using a 
charge count value of single-muon peak. 


(2) Scintillation counter: In order to have a certain 
position resolution of high energy particles in the air 
shower cores, the size of an inner detector unit is taken 
to be 80 cm x 50 cm, as shown in Fig. [3] One unit 
scintillation counter is divided into 20 pieces of 4 cm 
width. For getting a better uniformity of light output 
when a high energy particle hits different position of a 
scintillation counter unit, 40 WLSFs are installed paral- 
lelly in the center of each scintillator and connected to 
high-gain PMT (R4125) while 20 WLSFs installed in the 
center of each scintillator are connected to low-gain PMT 
(R5325). Such design guarantees geometrical uniformity 
of the detector response within 6%. 

(3) PMT: In order to record the electromagnetic 
showers in the energy range from several GeV to sev¬ 
eral 10 TeV, a wide dynamic range from 10 MIPs to 
10 6 MIPs of PMT is required. In addition, taking into 
account the importance of single-particle measurement 
in the system calibration, the dynamic range of PMT 
should be from 1 MIP to 10 6 MIPs. This is realized by 
adopting a high-gain PMT (HG-PMT, HAMAMATSU 
R4125) and a low-gain PMT (LG-PMT, HAMAMATSU 
R5325) that are responsible for the range of 1—3 x 10 3 
MIPs and 10 3 — 10® MIPs, respectively. 
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(4) WLSFs: To save the wave length shifting fibers, a 
short length (110 cm, half of the fiber length used in the 
old detectors 0 ) of fiber is used for which one end is 
connected with the PMT and another end is plated with 
aluminum for the reflection, as shown in Fig. [3] Since the 
reflection rate reaches 99% and only the gain of PMT is 
recorded from each detector, there is no influence on the 
experimental data. Considering the large temperature 
difference in Yangbajing within one year even one day, 
a heat-insulator layer is used inside each detector unit 
box. 

YAC-II array has been operating together with the 
Tibet-Ill array and the underground MD array since 
March, 2014, as shown in Fig. [T| The Tibet-Ill array 
is used to measure the shower size and the arrival di¬ 
rection of each air shower. Any four-fold coincidence of 
the FT detectors is used as the trigger condition for air- 
shower events. The air-shower direction can be estimated 
with an error smaller than 0.2° above 100 TeV [l|, and 
the primary energy resolution is estimated to be 12% at 
energies around 10 15 eV by our simulation [l3 |. 

Trigger rate of YAC-II is about 3.5 Hz with a dead¬ 
time rate of ~ 1%. If one YAC-II detector makes a 
trigger signal, all ADC data from all YAC-II units are 
recorded. Also the trigger signal is sent to DAQ sys¬ 
tem for Tibet-Ill array and MD array. ADC modules of 
YAC-II are calibrated every 4 hours. ADC pedestal val¬ 
ues are measured every 10 minutes. Each DAQ system 
has a GPS clock module independently. The matching 
between YAC, Tibet-Ill and MD data is made by taking 
a coincidence of GPS clocks, trigger tag to Tibet-Ill and 
MD array. The coincidence interval of GPS is shorter 
than 1 ms. The average time difference is 8.1 /zs ± 0.4 

fjiS. 


3 Detector Calibration 


3.1 The probe calibration 


In this hybrid experiment, the burst size is defined as 
the PMT output (charge) divided by that of the single¬ 
particle peak, which is determined by a probe calibration 
using cosmic rays, typically muons. For this purpose, a 
small scintillator 25 cm x 25 cm x 3.5 cm with a PMT 
(H1949) is put on the top of each detector during the 
maintenance period. This is called a probe detector, as 
shown in Fig. [4j Fig. [5] shows the charge distribution of 
single particle in a detector unit. The peak is defined as 
one MIP ( 1 MIP = 2.12 ± 0.01 pC ). Therefore, the 
determination of the burst size N b is calibrated by using 
a single-particle peak. 



Figure 4. The view of the probe calibration. 



Figure 5. The charge distribution of a single muon 
in the first detector unit. The peak called single¬ 
particle peak is defined as one MIP. The single¬ 
particle peak of this detector is 2.12 ± 0.01 pC. 


3.2 Linearity of PMTs 


For each PMT (R4125 and R5325) used in YAC-II, 
the linearity with respect to high voltage is better than 
1%, as shown in Fig. [G] The high-gain PMT (R4125) is 
used to measure light yield from 1 MIP to a few thousand 
MIPs, while low-gain PMT (R5325) cover the dynamic 
range from a few hundred to 10 6 MIPs. The dynamic 
range of PMTs is measured using LED light source and 
optical filters. In the test experiment we fixed the posi¬ 
tions of LED, filter and PMT. By using different filters, 
we can get lights of the different intensity. The LED is 
driven by TTL pulse with the width of 25 nsec. The 
details are described in 13]. 


3.3 Uniformity of the detector 


When the charged particles pass through the scintil¬ 
lator in different position, the light transmitted to the 
PMT by the WLSFs will vary in a little range which 
is called “position dependence” or “uniformity” of the 
detector. The uniformity of YAC-II detector has been 
measured using cosmic-ray single muons selected by a 


XXXXXX-4 


















Chinese Physics C Vol. XX, No. X (201X) XXXXXX 


triple coincidence as shown in Fig. [TJ As shown in the 
figure, the uniformity in the scintillator length 50 cm is 
better than 6%. 
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Figure 6. The linearity between the output (pC) 
and the high voltage supplied to high-gain PMT 
(R4125, Red) and low-gain PMT (R5325, Blue). 
The points show the experimental data. The red 
dotted line is a fitting line. 
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Figure 7. Top: Experiment set-up for single muon 
measurement. Bottom: the uniformity of YAC-II 
detector. 


3.4 Linearity of PMT plus scintillator 

The linearity and the saturation of the plastic scintil¬ 
lator and PMT used in the YAC detector have also been 


calibrated by using the accelerator beam of the BEPCII 
(Beijing Electron Positron Collider, IHEP, China). The 
accelerator-beam experiment shows a good linearity be¬ 
tween the incident particle flux and YAC-II output below 
5x10® MIPs and the saturation effect of the plastic scin¬ 
tillator satisfies YAC detector’s requirement. The details 
are described in fl5| . 

3.5 The correlation between high-gain PMT 
and low-gain PMT 

Since we can not measure a single muon using low- 
gain PMT (R5325) directly, we calibrate the low-gain 
PMT output using the the correlation between high-gain 
PMT (R4125) and low-gain PMT (R5325) in their over¬ 
lapping region. After the detectors running consistently 
for about 10 days, we get the gains of low-gain PMT 
(R5325) and high-gain PMT (R4125) of the same burst 
events which are recoded by ADC, as shown in Fig. [8] 
From Fig. [HI it is found that there is a fine linear rela¬ 
tionship between two PMTs data. We can then obtain 
the ratio (just the slope of the line in Fig. 0) between 
them. With a further calculation it may be possible to 
get N b only from output signals from low-gain PMT. 


Graph 



Output of HG-PMT(pC) 


Figure 8. The correlation between output signals 
from high-gain PMT (R4125) and low-gain PMT 
(R5325) obtained during 10 days operation of 
YAC-II. 


3.6 Consistency of the gains of different detec¬ 
tor units 

Setting a lower detection threshold (say, taking 30 
mV as the discrimination threshold) and “any 1” (at 
least any one unit of the 124 units should be fired) as 
the trigger condition, YAC-II was operated for 10 days 
to measure the W-spectrum of each detector unit. Some 
differences are found from the W-spectrum of each unit. 
By slightly adjusting the working voltage of PMTs, con¬ 
sistent spectra are obtained. In Fig. 0 it is shown that 
the difference of the W-spectrum between 124 YAC de¬ 
tectors is less than 20%. 
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Figure 9. The JV^-spectrum with the experimental 
data of YAC-II in 10 days. It is shown that the 
difference of the AVspectrum between 124 YAC 
detectors is less than 20%. 


4 Summary 

A new air-shower core detector (YAC-II) has been set 
up at Yangbajing and started data taking from March, 
2014. Each detector of YAC-II consists of 0.4 m 2 plastic 
scintillator, equipped with a high-gain PMT and a low- 
gain PMT which is put under 3.5 cm lead plate. The 
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